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ABSTRACT 

We report a three year, 18 epoch, VLBI monitoring study of H2O masers in the sub-parsec, 
warped, accretion disk within the NGC 4258 AGN. Our immediate goals are to trace the geometry 
of the underlying disk, track rotation via measurement of proper motion, and ascertain the radii of 
masers for which centripetal acceleration may be measured separately. The monitoring includes 
~ 4x as many epochs, '-^ 3x denser sampling, and tighter control over sources of systematic 
error than earlier VLBI investigations. Coverage of a ~ 2400 kms~^ bandwidth has also enabled 
mapping of molecular material 30% closer to the black hole than accomplished previously, 
which will strengthen geometric and dynamical disk models. Through repeated observation we 
have also measured for the first time a 5/ias (Icr) thickness of the maser medium. Assuming 
this corresponds to the thickness of the accretion disk, hydrostatic equilibrium requires a disk 
plane temperature of « 600 K. Our long-term goal is a geometric distance to NGC 4258 that is 
accurate to ~ 3%, a ~ 2x improvement over the current best estimate. A geometric estimate of 
distance can be compared to distances obtained from analysis of Cepheid light curves, with the 
intent to recalibrate the extragalactic distance scale with reduced systematic uncertainties. This 
is the first paper in a series. We present here VLBI observations, data reduction, and temporal 
and spatial characteristics of the maser emission. Later papers will report estimation of orbital 
acceleration and proper motion, modeling of disk 3-D geometry and dynamics, and estimation of 
a " maser distance." Estimation of a "Cepheid distance" is presented in a parallel paper series. 

Subject headings: galaxies: individual (NGC, 4258) — galaxies: kinematics and dynamics — galaxies: 
nuclei — galaxies: active — masers — techniques: interferometric 



1. Introduction 

The extragalactic distance scale (EDS) is un- 
derpinned by measurements of "standard candles" 
in the Large Magella nic Cloud (LMC ) and other 
nearby galaxies ( Freedman et al. .I I2OOII ). The long- 
term importance of a high-accuracy EDS calibra- 
tion lies in potential precision estimation of cosmo- 
logical parameters, su ch a,s the equation of state 



for dark energy (e.g., |HuJ 120051). One means by 



which the EDS calibration might be improved is 
by anchoring "standard candle" calibration to dis- 
tances estimated by geometric techniques. To- 
wards this end, NGC 4258 — for which a geometric 
distance has been estimated from modelling of the 
3-D geometry and dynamics of the accretion disk 



surro unding the central engine (jHerrnstein et al 
19991) — could either replace or be used in conjunc- 
tion with the LMC as anchor. 

In NGC 4258, accretion disk material emits 
H2O maser emission in the 616 — 523 transition at 
22.235 GHz, which was fir st reported close to th e 
systemic velocity iVsys) bv lClaussen et al. ( 1984 ). 
A key discovery was high-velo city Doppler compo- 
nents at ~ ^^^^±1000 km s~^ ( Nakai. Inoue. fc Mivoshil 



r sys 



472 ± 4 



I993I ) in NGC 4258, where 
km s~^, adopting the Local Standard of Rest 
(LSR) and radio definition of Doppler shift 



(jCecil. Wilson, fc Tullvl Il992h . An early Very 



Long Baseline Interferometry (VLBI) study l| Greenhill et al 



^The VLB A is a facility of the National Radio Astronomy 
Observatory, which is operated by Associated Universities 
1 Inc., under cooperative agreement with the National Sci- 
ence Foundation. 



1995al ) suggested that maser emission might orig- 
inate from a rotating disk, as did a parallel in- 
vestigation of secular velocity drifts for the peaks 
of spectral features, i.e., centripetal accelerations 
(Haschick, Baan, & Peng 1994; Greenhill et al. 
1995b; Nakai et al. 1995; see also Watson & Wallin 
1994). Low- velocity Doppler components (those 
near Vsys) drifted by ~ -1-10 km s~^ yr~^ and high- 
velocity components drifted by < ±lkms~^ yr~^, 
thus localizing the former to the near side of 
the disk and the latter near the mid-line (i.e., 
the diameter perpendicular to the line of sight). 
The sub-parsec disk model was confirmed by 
Very Long Baseline Array (VLBA)0 observations, 
which traced, most importantly, the disk warp 
and the rotation curve, the latter being Keplerian 
to better than 1% (Miyoshi et al. 1995; see also 
Herrnstein et al. 2005). 

For a Keplerian system, accurate first-order 
estimates of geometric distance may be ob- 
tained relatively simply from the geometry of the 
warped disk and measurements of acceleration or 
proper motion for low-velocity emission. To date, 
NGC 4258 is the only galaxy for wh ich this has 
been done. Herrnstein et al. ( 1999( ) obtained a 
distance to NGC 4258 of 7.2 ± 0.3 {random) ± 
0.4 (systematic) Mpc, where the systematic term 
is due mostly to a thus far weak constraint on 
orbital eccentricity (< 0.1). Our major aim is 
to reduce both rando m and systematic error. 
Herrnstein et al. ( 1999f ) used four epochs of VLB A 
data spanning 3 years to estimate distance. In 
contrast, we have 18 epochs over a subsequent 3.4 
year interval. This facilitates more detailed mod- 
eling of disk substructure, formal incorporation of 
orbital eccentricity and periapsis angle as param- 
eters, more definitive tracking of disk rotation, 
and reduction of randon i and systematic e r ror by 

(|2005h 



Humphreys et al 



a factor of at least two. 
discuss early progress. 

The current distance scale, as characterized by 
the Hubble constant (Hp) , may be uncertain b y 
perhaps as little as 10% ( Freedman et al. 200l[ ). 
The most robust estimate depends on measured 
mean magnitudes of ^ 800 Cepheid variable stars 
in 31 galaxies within 30 Mpc, and comparison with 
the period luminosity (PL) relation for Cepheids 
in the Large Magellanic Cloud (LMC). The LMC 
establishes a zero point for calibration of Cepheid 
light curves in other galaxies. Common secondary 



distance indicators that are identifiable in the 
Hubble Flow (e.g., TuUy Fisher relation, type la 
SNe), are calibrated using distances to Cepheid 
host galaxies and are used in turn to estimate 
Ho . High-precision estimates of Ho are also avail- 
able from analysis of Cosmic Mic rowave Back- 



groun d (CMB) fluctuations (e.g., ISpergel et al 



20061 ). However, degeneracies in analysis of CMB 



data do not allow independent determination of 
both Ho and the equation of state for dark energy, 
w. Note that Spergel et al. ( 20061 ) assumed a flat 
universe with w = —1 and obtained a value of Ho 



in good agreement with iFreedman et al.l (|200ll ). 

The most substantial problem with the current 
Cepheid-based calibration of the EDS is depen- 
dence on the LMC as zero point anchor for PL re- 
lations. Contributing factors include uncertainty 
in distance to the LMC, poorly understood inter- 
nal galactic structure, and controversy over the 
impact of the difference in stellar mctallicity be- 
tween the LMC and HST Cepheid-sample galax- 
ies. 

The adopted unc e rtaint y in the LMC distance 
by iFreedman et al. (2001) is '~ 5%. However, 
LMC distance estimates obtained from different 
calibrations and analyses over the last decade 
cover a ±23% range (±0.50 mag when expressed 
in terms of distance moduli). This is well beyond 
unc ertainties for ind ividual measurements (Figure 
8 of lBenedict et al.1[2002a) and likely indicative of 
unresolved systematic errors. Although uncertain- 
ties ^ 10% have been obtained from the distribu- 
tions of distance estimates for different subsam- 



ples of measurements (e.g., IFreedman et al.l 12001 



Benedict et al.l l2002bl ) . in the face of apparently 



unresolved systematics, the averages and moments 
may be statistically problematic. 

The uncertain structure of the LMC (e.g., 
disk eccentricity, out-of-plane stellar component), 
complicates interpretation of distance in dicators 
van der Mareill200ll : lNikolaev et all2004 and ref- 



erences therein). For instance, robust estimates 
of distance have been obtained for eclipsing bina- 
ries, with individual uncertainties of 2 - 3% (e.g.. 



Fitzpatrick et al.l I2002L and references therein). 



However, for the four systems well studied thus 
far, the range of i mplied distances for th e LMC 
center covers 14% ([Fitzpatrick et al.ll2003l) . 

The impact of metallicity could be as large at 
±10% in distance. Limited theoretical and cmpiri- 



2 



cal understanding has led to long term controversy 
over the magnitude and sign of the effe ct of metal 



Groenewegen et al.i2004 : Sakai et al.l2004 Romaniello et 



conte n t on Ce pheid brightness ( Sa sselqv et al 
1997t iKochmekl Il997l: iKennicutt et all Il998l: 



maser, i.e., three 200 km s ^-wide lin e complexes 
distributed over ^ 2000 kms~^ (e.g.. iBragg et al 



2OOOI ): (3) tight control over sources of systematic 



20061 ). Attempts to make definitive measure- 
ments have been complicated by the magnitude 
of observable metallicity gradients in well studied 
galaxies, the sizes of Cepheid samples, and the ef- 
fects of reddening and crowding in low-metallicity 
fields at smaller galactocentric radii. Consensus 
theoretical models have not emerged because of 
computational difficulty, with some doubt as to 
whether or not metal content alone needs to be 



considered (jFiorentino et al.l 120021 ). 

Precision estimation of cosmological par ame- 
ters is best achieved by multiple approaches. 
([20051) argues that percent level accuracy in Ho, 
combined with measurements of CMB fluctua- 
tions, would enable robust constraint on w in- 
dependent of those provided by measurements 
of, e.g., weak lensin g and large scale structure. 



Spergel et al. 



( 20061 ) also demonstrate the im- 
pact on estimates of curvature by independent 
constraint on Hq. A new "maser distance" for 
NGC4258 is a first step, with the immediate goal 
of ^ 3%. This paper is the first in a series. Here 
we present the VLBI data and briefly discuss gen- 
eral trends such as time evolution of position- 
velocity traces. Later work will present time-series 
maser spectroscopy and modeling of the combined 
datascts to o btain an improved geometric distance 
(cf . |Hcrrnstcin_et_alj|T999) . Parallel analysis of re- 
cent Cepheid photometry obtained with the Hub- 



ble Sp ace Telescope, beginning with lMacri et al. 
(I2OO7I ). will provide an improved "Cepheid dis- 



yuuxlj, will provide an imp : 
tance" (cf. iNewman et al. I 



20011 ). Ultimately, 



the two programs will culminate in comparison 
of maser and Cepheid distances and re-assessment 
of the EDS and implications for cosmology. 

2. Observations and Data Reduction 
2.1. The Observations 

The design of the observing program was in- 
tended to enable (1) unambiguous tracking of 
evolution in the source spectrum and disk struc- 
ture, through close spacing of epochs; (2) monitor- 
ing to detect possible transient emission outside 
the historically recognized velocity range of the 



er ror; and (4) relatively uniform sensitivity from 
adit och to epoch. 

We observed NGC4258 18 times between 1997 
March 6 and 2000 August 12, using the VLBA, the 
Very Large Array (VLA), and the Effelsberg 100- 
m telescope (EFLS) of the Max-Planck-Institut 
fiir Radioastronomic (Table 1). Twelve "medium- 
sensitivity" epochs involved the VLBA alone; the 
remaining six "high-sensitivity" VLBI epochs were 
interleaved and involved the VLBA, augmented by 
the phased VLA and EFLS. 

In medium-sensitivity epochs, we recorded 
eight 16 MHz (« 216 km s^^) basebands with 
1-bit sampling. We overlapped seven bands by 
15%, resulting in observing bandpasses of ~ 362 
to 1676 km s~^ (covering low- velocity and red- 
shifted high- velocity emission) and ~ —706 to 
608kms~^ (covering low- velocity and blue-shifted 
high- velocity emission), depending upon the epoch 
(Table 2). We tuned the eighth band to enable 
dual polarization observations of the emission 
peak (spanning the velocity range 362-578 km 
for red-shifted epochs and 392-608 km s~^ for 
blue-shifted epochs), thus increasing the signal-to- 
noise ratio available for tracking atmospheric and 
instrumental response fluctuations during calibra- 
tion. Because the ~ 2000 km s^^ bandwidth of 
the source was difficult to accommodate instan- 
taneously with VLBI data acquisition systems, 
we alternated month-to-month between observ- 
ing two offset ^ 1500 kms^^ bands. One covered 
low-velocity and red-shifted high-velocity maser 
emission, and the other covered low-velocity and 
blue-shifted high-velocity emission. Combining 
data from adjacent months, we effectively ob- 
served ~ 2400 km s^^, centered approximately on 
the systemic velocity. 

In high-sensitivity epochs, we recorded eight 8 
MHz (« 108kms~^) basebands with 2-bit sam- 
pling, alternating baseband tuning over the course 
of each track to cover the three emission com- 
plexes. The narrower instantaneous bandwidth 
accommodated the limited IF bandwidth of the 
VLA (2x50 MHz) . In the first five high-sensitivity 
epochs, we recorded both upper and lower side- 
bands to limit the number of baseband mixers re- 
quired and thereby accommodate a reduced com- 
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plement at Effelsberg. To avoid band gaps at 
mixer frequencies, we overlapped the sidebands of 
adjacent mixers and obtained « 300kms~^ cov- 
erage of low- velocity, red-shifted, and blue-shifted 
emission complexes. In the last high-sensitivity 
epoch, we took advantage of a new complement 
of mixers at Effelsberg, recording only upper side- 
bands with 20% overlap and 365kms~^ coverage 
of each emission line complex. 

2.2. Strengths of the Current Study 

Time Sampling — Spacing between high-sensitivity 
epochs, ^ 8 months, was comparable to spac- 
i ng between epochs in an earlier VLBI study 
(jHerrnstein et al.l bOOSl their Table 1), but the 
1-2 month time interval between most medium- 
sensitivity epochs was comparable to the typi- 
cal known t imescale for spectral variability (e.g., 
Bragg et al. lHooo). A minimum 8 day separation 
between epochs also enabled testing for short-term 
source structure variability. 

Bandwidth — The high-sensitivity epochs (as in 
all previous VLB I studies) focused on the velocity 
ranges in which Nakai. Inoue. fc Mivoshil ( 19931 ) 
identified emission. In contrast, medium sensi- 
tivity observations monitored all velocities from 
Vlsr ~ —706 to 1676 km s~^, achieving a detection 
threshold « 2x below that of earlier broadband 
(—1200 to 2650 km s~^) singl e-dish studies by 
Nakai. Inoue. fc Mivoshi ( 19931 ) and Nakai et al 



median for medium-sensitivity epochs (Table 1). 
Earlier studies of low- and high-velocity emission 
together achieved s ensitivities of 3.0 - 6.5mJy with 
a 5.6 mJy median ( Herrnstein et al.ll2005l ). 



(|1995( ). scaled to the same channel spacing. Emis- 
sion at "new" velocities was sought because it 
would presumably arise in previously unmapped 
portions of the accretion disk. For each epoch, we 
recorded eight basebands (which were correlated 
in two passes). Each baseband was divided into 
512 spectral channels. 

Quality and Consistency — The same observ- 
ing setup and data reduction path were used for 
each medium-sensitivity epoch. A separate com- 
mon setup and reduction path was used for high- 
sensitivity epochs. Deliberate repetition (e.g., u,v- 
coverage, choice and timing of calibrators, base- 
band tuning) limited systematic errors in differ- 
ential measurements of dynamical quantities (i.e., 
proper motion, acceleration). 

Sensitivity — Noise levels (lu) in synthesis images 
were 2.3 - 4.7 mJy with a 3.0 mJy median for high- 
sensitivity epochs and 3.6 - 5.8 mJy with a 4.3 mJy 



2.3. Calibration and Imaging 

Calibration and synthesis imaging were per- 
formed using the Astronomical Image Processing 
System (AIPS) software package with largely stan- 
dard techniques for spectral-line VLBJ^. We dis- 
cuss notable elements of the processing below. 

2.3.1. Array and Sky Geometry 

We sought to ensure consistency in the an- 
tenna coordinates used for each correlation over 
the 3 years of monitoring. This was achieved by 
implementing (post-processing) the station posi- 
tions and (tectonic) velocities estimated by the 
United States Naval Observatory (USNO) in so- 
lution N9810 (Table 3). Because the VLA does 
not have dual-frequency geodetic VLBI receivers 
and is absent from solution N9810, we adopted 
a position estimated by the NRAO using single- 
frequency data for epoch 2000.9 and the average 
velocity for the nearby Pietown and Los Alamos 
stations (C. Walker, private communication). Af- 
ter correction, station positions were accurate to 
w 1 mm for the VLBA and EELS and w 3 cm for 
the VLA. 

To ensure consistency in the maser posi- 
tion and error budget from epoch to epoch, we 
adopted a maser position obtained from analysis 

of early VLBI data: a(J2000) = 12^18^57 f 5046 ± 

Of 00 03: (5(J2000) = 47°18'14:'303±0:'003 (iHerrnstein et al 
20051 ). This was used in correlation of all but 



the first three epochs, for which we applied a 
frequency-dependen t shift to the measured maser 
spot positions, after Herrnstein et al. ( 20051 ). The 
shift corresponded to an error of —0 50054 and 
0''033 in right ascension and decHnation, respec- 
tively, in the a priori position. A priori calibrator 
positions were accurate to < 1 mas, and no cor- 
rections were made. A priori earth orientation 
parameters were compared to "final" values pub- 
lished by the USNO and found to be accuarate to 
< 1 mas (pole position in X and Y) and < 0.1 ms 



^AIPS cookbook, chapter 

http: / /www, aoc.nrao.edu/aips/cook. html I 
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(UTl-UTC) for each epoch. We discuss the im- 
pact of these uncertainties and hmits in §2.5 (see 
also Table 4). 

2.3.2. Amplitude Calibration 

To calibrate VLBA fringe amplitudes, we ap- 
plied observatory-standard antenna gain curves 
and time-series system temperature measure- 
ments. An w 40% depression in the system tem- 
peratures of the red-shifted high-velocity base- 
bands in high-sensitivity epochs rendered some of 
these calibrations suspect. Such large changes in 
system temperature is not likely and is probably 
due to a measurement problem. To compensate, 
we computed a system temperature time-series 
characteristic of the several low-velocity base- 
bands and applied it to all basebands. 

Accurate system temperature time-series were 
typically not available for the VLA and EFLS. 
For these two stations, we calibrated amplitudes 
through a comparison of fringe amplitudes on 
4C39.25 among VLBA stations (after calibration) 
and fringe amplitudes for baselines including the 
VLA or EFLS. Time sampling of the resulting cal- 
ibration was relatively coarse, but we anticipated 
later self-calibration would enable refinement. 

We note that in the three cases for which we 
made nonstandard amplitude corrections (i.e., 
VLA, EFLS, and the red-shifted high- velocity 
VLBA basebands in high-sensitivity epochs) am- 
plitude errors for some basebands could survive 
self-calibration, because gain corrections are cal- 
culated for emission in just one (low-velocity) 
baseband. Though a matter of concern, ampli- 
tude errors on the order of a few tens of percent 
on a small fraction of baselines are anticipated to 
have a relatively small impact on astrometry for 
point-like sources, i.e., maser spots. 

We expect the overall flux scale for a given 
epoch to be accurate to ~ 30% for medium- 
sensitivity epochs and ~ 50% for high-sensitivity 
epochs. However, when comparing the overall flux 
scale between or among epochs, we expect an accu- 
racy of « 15% for medium-sensitivity epochs and 
w 25% for high-sensitivity epochs. 

2.3.3. Global Fringe fitting 

The most demanding element of the global 
fringe fitting process was the accurate measure- 



ment of the slowly time variable electronic phase 
offsets among basebands, using observations of 
quasars made approximately every hour for that 
purpose. Specifically, the wide spacing in fre- 
quency among bands increased susceptibility to 
(nonphysical) phase wraps entering the interpo- 
lation of phase offsets between adjacent calibrator 
scans. Because the phase wraps were readily de- 
tected, straight-forward editing of single-band de- 
lay and rate data enabled removal and estimation 
of offsets with accuracies on the order of 10° or 
3% of a fringe on the sky. 

2.3.4- Self- calibration 

We applied self-calibration solutions to the data 
in each baseband to correct for fluctuations in 
fringe amplitude and phase. These fluctuations 
arose from instabilities in the atmosphere and in 
the frequency standards at the stations and from 
antenna gain changes. For each epoch, we identi- 
fied a relatively unblended Doppler component (^ 
1 Jy) to provide the reference signal. Images of the 
reference component achieved dynamic ranges of 
40-300, defined with respect to the 5a noise level 
or strongest artifact. 

We observed unusually large artifacts in images 
of high- velocity red-shifted emission for a few high- 
sensitivity epochs. Dynamic ranges were reduced 
by a factor of ^ 4. We were unable to identify a 
cause and adopted an empirical strategy to cor- 
rect the problem. For these epochs, we performed 
a second self-calibration on the strongest observed 
red-shifted Doppler component. Though typically 
only a few tenths of a Jy, this secondary reference 
was detectable because the initial self-calibration 
extended the phase coherence of the data. We 
applied the additional, slowly varying amplitude 
and phase corrections to the data in all red-shifted 
basebands. To recover relative position informa- 
tion (with respect to the primary reference emis- 
sion), we applied the angular offset between refer- 
ence features, measured prior to the second self- 
calibration. 

Secondary self-calibration of the blue-shifted 
high-velocity emission was not possible because 
the peak flux density was too low (< 100 mJy). 
We tested two secondary calibrations for the "blue 
data". The first was achieved by applying the 
"red data" self-calibration corrections to the "blue 
data" and the second by applying the complex 
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conjugate of the "red data" corrections to the 
"blue data" . In the first case we implicitly mod- 
eled the unexplained errors as phase-only, while in 
the second case as delay-only. Neither approach 
improved the signal-to-noise in images of blue- 
shifted emission, and was therefore not adopted. 
Some losses in imaging blue-shifted emission re- 
main possible. 

2.3.5. Imaging and Deconvolution 

We imaged and CLEANed a 25.6 x 25.6 mas 
field (in medium-sensitivity epochs) or a 30.72 x 
30.72 mas field (in high-sensitivity epochs) over 
the inner 96% of each bandpass. The known angu- 
lar extent of the maser source is ~ 1 x 16 mas. Syn- 
thesized beam sizes and typical noise levels for in- 
dividual spectral channels are given in Table 1. We 
used a hybrid weighting scheme midway between 
uniform and natural weighting (ROBUST=0 in 
the AIPS task IMAGR), in order to balance the 
desire for a small beam measured at the half-power 
points and low rms noise. In five of the six high- 
sensitivity epochs, we extended imaging outside 
the inner 96% of the bandpasses, because a blue- 
shifted Doppler component had been found in the 
medium-sensitivity epochs. The anticipated emis- 
sion was detected in one epoch (BM081A). 

2.4. Identification of Maser Spots and 
Continuum Emission 

We fit a 2-D Gaussian brightness distribution to 
all peaks in channel maps and required that they 
satisfy the following requirements, in order to be 
judged bona fide maser "spots." 

• Emission centroids lie within one beamwidth 
over at least three adjacent channels at or 
above the 6a level. 

• Fitted peak intensities exceed the intensity 
of imaging artifacts, i.e., intensities must ex- 
ceed the local peak intensity divided by the 
dynamic range. We note that this dynamic 
range cut was only relevant to channels in 
high-sensitivity epochs with strong emission 
features and candidate secondary features. 

• Half-power major and minor axes are both 
smaller than twice the major and minor axes 
of the convolving beam. In cases where 



breadth appeared consistent with blending 
of individual maser spots, we attempted to 
fit multi-component Gaussian models. 

In each epoch, there are on the order of 
10^ independent pixels among the image cubes. 
The probability that one or more will exceed 
6(7, in the absence of a signal, is « 10% (e.g., 
iThompson. Moran. &: Swe nson 200l Fig. 9.5). 
By requiring three adjacent channels to satisfy 
this criterion, with a constraint on position, the 
chances of false detection are on the order of 10~^ 
(where we have accounted for correlation among 
channels consistent with the FX architecture and 
time- windowing of data in the VLBA correlator) . 

Where we had a priori information as to the 
location in velocity and position of a maser, we re- 
laxed the 6(7 selection criterion. There were three 
instances. 

• High-velocity features (> 3.5(7) were ac- 
cepted if emission at corresponding po- 
sitions and velocities had been detected 
above the formal detection threshold in 
other epochs. The probability of false de- 
tection of a 3.5(7 emission peak over one 
independent point (one beam) is < 10~^. 
Blue-shifted Doppler components falling into 
this category were found at —285 km 
(BM081A), -367 km s'^ (BM112G), -373 
km s-i (BM081B and BM112C), and -515 
km s-i (BM112C, BM112E, and BM112H). 
Red-shifted Doppler components falling into 
this category were found at 1566 km s"^ 
(BM112F and BM112M). 

• High- velocity blue-shifted emission (> 2a) 
was accepted in spectral channels adja- 
cent to already identified blue-shifted fea- 
tures, if in the same epoch and at the 
same angular position to within a beam 
width. Blue-shifted emission was identi- 
fied in this way at -285 km 8"^ (BM112C, 
BM112E, BM112G, BM112L, BM112N, 
and BM112P), -367 km s^^ (BM112C), 
-373 km s-i (BM056C, BM112E, BM112G, 
BM112H, and BM112L), -435 km s^^ 
(BM056C and BM081A), -441 km s^^ 
(BM056C), and -515 km s'^ (BM112G). 
The motivation for this extension of our de- 
tection limit was better establishment of line 
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shape for weak features. It was not used to 
increase the number of features. 

High-velocity featu res (> 3.5a) outside the 
nomi nal range fi.e.. lNakai. Inoue. Sz Mivoshi 



19931 ) were accepted, if single-dish study sep- 



arately identified the Doppler components. 
In this way, we detected th e 1647 km s~^ 
compo nent first reported by iModi azetaL 
(|2005l ) for epoch 2003.75. Extrapolat- 



ing from the positions of other red-shifted 
Doppler components and adopting a Kep- 
lerian rotation curve for the "maser disk," 
we confined the anticipated position of the 
1647 km s~^ emission to within 1 beam 
X 3 channels. Emission was detected in 
BM112M at 5(7 in the line-center channel 
and 4a in adjacent channels. Marginal de- 
tections were also found in two adjacent 
epochs, i.e., emission at 2.4 and 3.0(7 in two 
contiguous spectral channels for BM112K 
and four 2.8(T — 3.3(T points in adjacent spec- 
tral channels for BM1120. 

In addition to identifying maser emission, we 
searched for continuum emission at all epochs 
within 12.5 mas of the low- velocity reference emis- 
sion. Continuum u,v-data were constructed by fre- 
quency averaging each baseband prior to imaging. 
In the case of BG107 and all medium-sensitivity 
epochs, we discarded the outer 2% of each of the 
red- or blue-shifted basebands and obtained an ef- 
fective average of 1960 channels (30.6 and 61.3 
MHz for BG107 and medium-sensitivity epochs, 
respectively). In the other high-sensitivity epochs, 
we discarded the outer 2% of each of three red- or 
blue-shifted basebands and averaged 1470 chan- 
nels (23.0 MHz). Fewer channels were retained in 
these epochs, because of extensive overlap of side- 
bands for adjacent baseband mixers. 

Finally, we performed a separate search for 
maser emission that might be associated with 
jet activity in NGC4258, well downstream from 
the central engine. In the most sensitive of our 
high-sensitivity observations (BM112H), we box- 
car averaged the u,v-da,ta, over four channels to 
improve sensitivity and imaged a 30.72 x 491.52 
mas (512 x 8192 pixels, E-WxN-S) strip centered 
on the reference position. No emission peaks ex- 
ceeded a 6a detection threshold (13mJy). 



2.5. The Error Budget 

For each Gaussian fit, the random measurement 
error in position was estimated as being the larger 
of the fitting error or the beam full-width at half- 
maximum (FWHM) divided by twice the peak 
signal-to-noise ratio {S/N) . We estimated the sys- 
tematic position error from two components: (1) 
errors that scale with difference in sky frequency 
between a given channel and the reference (i.e., 
v — t're/)) and (2) errors that scale with the dif- 
ference in baseband frequencies of a given channel 
and the reference (i.e., [v — Vo] — [vref — i^ref^]: 
where Vo and t're/o are the band-center frequen- 
cies). The maximum difference in sky frequencies 
is « 75 MHz for NGC4258; the maximum differ- 
ence in baseband frequencies is < 8 MHz. Er- 
rors in the absolute position of the reference and 
calibrators, or errors in zenith propagation delay 
and baseline coordinates give rise to the first type 
of systematic error. Uncertainty in station clocks 
give rise to the second type. 

A consistency check of the high-accuracy (^ 3 



mas) m aser reference position of lThompson. Moran. fc Swenson 
(l200ll) can be obtained from our continuum imag- 
ing of BM081B. In this epoch, we imaged with a 
less accurate (~ 60 mas) a priori maser reference 
position, which led to a substantial difference be- 
tween red- and blue-shifted position offsets, due 
to the large frequency difference between red- and 
blue-shifted continuum "bands" (see Table 4). If 
we assume that this misalignment of red- and 
blue-shifted continuum is entirely due to an inac- 
curate a priori maser reference position, we ob- 
tain a new reference position that is offset from 
Thompson. Moran. fc SwensonI (|200l[) by OfOOl in 
right ascension and 0"01 in declination, well within 
our 1.5a combined errors (ours plus Herrnstcin's, 
added in quadrature). We note that all channel 
and continuum peak positions in this epoch were 
corrected, post-imaging, for the lower-accuracy a 
priori maser reference position. 

Table 4 gives the magnitudes of random and 
systematic position error. Random error domi- 
nates the error budget for weak masers, whether 
they be at low- or high- velocities. For instance, 
blue-shifted masers are dominated by random 
error. Systematic error dominates red-shifted 
masers stronger than w O.lJy and low- velocity 
masers stronger than « 0.5 Jy. The largest con- 
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tribution to systematic error arises from uncer- 
tainty in the absolute position of the reference 
emission, causing « 10 /xas error across 75 MHz, 
the typical frequency separation between low and 
high- velocity masers. For a given Doppler compo- 
nent, the error repeats from epoch to epoch, with 
second-order impact (because reference emission 
velocities differ by < 4 MHz among epochs) . How- 
ever, a 10 /xas error is comparable to upper limits 
on disk thickness an d the hydrostatic scale height 



tain a relativistic LSR velocity, 

t^rtop + Awlsr 



(|Moran et al.l Il995f ) , with implications discussed 



in a later paper. 

2.6. Relativistic Velocity Assignments 

Velocities in the LSR reference frame, assum- 
ing the radio definition of Doppler shift, were as- 
signed to each channel. The choices of reference 
frame and definition wer e made for cons istency 
with previous work (e.g.. iHerrnstein et al.l [20051. 
However, the first-order (radio definition) and full 
relativistic treatment of velocity differ by as much 
as 4kms~^ for the most red-shifted emission. Be- 
cause this is comparable in magnitude to impor- 
tant physical effects (e.g., distortion of the rota- 
tion curve by disk warping), we have adopted rel- 
ativistic velocities for purposes of disk modeling. 
Both first-order and relativistic velocities are pre- 
sented with the fitted maser spot parameters in 
Table 5. 

In order to convert classical to relativistic ve- 
locities, we first reconstructed the actual sky 
frequency of a feature observed in the earth's 
topocentric frame at 12'' UT for the date of ob- 
servation. This involved removing the LSR cor- 
rection to obtain the classical topocentric velocity 
and converting to the observed frequency accord- 
ing to the classical radio Doppler formula 



1^0 



= 1 - 



^ctop 



where Vohs is the observed sky frequency, vq the 
frequency in the emitter's rest frame, t^ctop is the 
classical velocity in the topocentric frame, and c 
is the speed of light. We then converted to a rela- 
tivistic topocentric velocity, 

Vrtop _ (1 - x'^) 
c ~ (1 + 2:2) ' 

Finally, we reintroduced the LSR correction to ob- 



WrLSR 



1 + (wrtop Ai;lsr)/c2 



After velocity correction, all positions were re- 
referenced to fiLSR = 510.0 km s^^, corresponding 
to the classical LSR velocity WcLSR = 509.6 km 
s^^ We note that in all previo us NGC4258 stud- 
ies (e.g., Herrnstein et al.l ll999V the quoted veloc- 
ities have been classically defined and conversions 
performed during accretion disk modeling made 
the analyses relativistically correct. Those con- 
versions included a general relativistic correction, 
which for a black hole of mass 3.8 x IO^Mq was 
4kms~^ at the disk inner radius and 2 km s~^ 
at the outer radius. We have not applied a gen- 
eral relativistic correction to the data presented 
here because the correction is model dependent, 
and its discussion is reserved for a later paper. 

3. Results 

3.1. Data Products 

Maser Spectra — Spectra of low- and high- 
velocity masers are shown in Figures 1-4. These 
spectra were constructed from fitted maser spot 
flux densities. An average spectrum for all epochs 
is shown in Figure 5. Channels without detectable 
emission are not shown. 

Maser Maps — We present the full set of Gaus- 
sian fits for detected maser spots in Table 5. Note 
that spots that lie in the overlap between adja- 
cent basebands appear twice. Quoted position 
uncertainties are fitted random error. The true 
random error is estimated to be the maximum 
of the fitted random error and the random er- 
ror computed according to Table 4. The sky dis- 
tribution and related statistics of the masers are 
shown in Figures 6-9; the position along the ma- 
jor axis of the disk versus LSR velocity are shown 
in Figures 10-13. We define the impact parame- 
ter to be 6 = ±{Ax^ + Ay2)i/2^ where Aa; and 
Ay are the east-west and north-south offsets from 
the emission at 510.0 km s""'^ (relativistic), respec- 
tively, and the sign of b is taken from the sign of 
Aa; (which is feasible because the disk is elongated 
chiefly east- west). The formal uncertainty in im- 
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pact parameter (ab) is given by 

(Aa; a^)^ + (Ay ay)^ 



1/2 



where ax and ay are the random components of 
uncertainty in east and north offsets, respectively. 

Continuum Maps — We detected continuum emis- 
sion > 3(7 in fifteen epochs and averaged decon- 
volved images (with uniform weighting) to obtain 
a peak intensity of l.f mJy {23a) (Figure f4). 
Nondetection for the remaining three epochs may 
be attr ibuted to known tim e variability of the 
source ([Herrnstein et al.lll998l ). In Table 6, we list 
(angle integrated) flux densities and centroid posi- 
tions. Figure 15 shows a plot of flux density versus 
time. Our measurements are consistent with the 
hypothesis that the emission traces a jet extend- 
ing outward from the centr al engine, along the spi n 



axis of the accretion disk ( Herrnstein et al. I ll998h . 



3.2. Position Shifts 

Visual comparison shows that the position- 
velocity loci for low- and red-shifted high-velocity 
emission are offset at a few epochs by up to 
~ 30/ias from a weighted average locus computed 
from the medium-sensitivity epochs. The shifts 
do not appear to be frequency dependent. We 
attribute them to residual phase offsets between 
basebands correlated in separate passes (as were 
low and high- velocity bands), uncertainty in the 
post self-calibration registration that was required 
for particular high-sensitivity epochs, and unmod- 
eled changes in the structure of the reference emis- 
sion from epoch to epoch. 

We computed corrections for this effect as fol- 
lows. First, we constructed variance- weighted av- 
erage X and y positions for each channel from all 
medium-sensitivity epochs. (Weighting depended 
upon random measurement error alone.) Next, 
we downweighted outliers according to the em- 
pirical scheme below and computed new variance- 
weighted averages. To downweight outlying sam- 
ples, we multiplied measurement errors by an em- 
pirical factor that increased with offset: exp(^)^, 
where A is the offset in a; or y from the weighted 
average and a is the unmodified measurement un- 
certainty. For instance, the measurement error for 
points that deviated by 6a is increased by a factor 



of 10, whereas the measurement error for points 
that deviated by la is increased by a factor of 1.06. 
Finally, we obtained offset corrections for individ- 
ual epochs (both high- and medium-sensitivity) by 
minimizing the sum of the squares of the offsets 
from the variance-weighted averages above. Fig- 
ure 16 shows the observed and corrected offsets for 
red-shifted emission in observation BM112H, an 
epoch with one of the largest shifts. Table 7 quan- 
tifies the computed shifts for each epoch, which we 
applied to the data. However, we note that Table 
5 contains uncorrected maser spot positions. In 
addition, no corrections have been computed for 
blue-shifted emission because it is too sparse to 
enable comparisons. 

4. The Sky Distribution of Maser Emis- 
sion 

Maser emission in NGC 4258 traces a thin. 



slight l y warped, nearly edge-q n disk (iMivoshi et al 
19951: (Herrnstein et al.l Il998l l2005f) . The new 
VLBI data (Figure 6) trace this structure more 
fully, due to both emission time variability (i.e., 
development of new Doppler components) and 
broader observing bandwidth (e.g., velocities 
above 1500 km s^^). The flaring (in north off- 
set range) of low-velocity maser spots at the east 
offset edges of the plot is an artifact and is due 
to the fact that the weaker spots, which have 
lower positional accuracy, are preferentially lo- 
cated away from the middle of the low-velocity 
spectrum (Figure 1). Crowding of plot symbols 
also obscures the relative breadths of the distribu- 
tions of maser spots with different flux densities. 
To assess the statistics of the distribution of maser 
spots, we have computed perpendicular offsets of 
low-velocity masers from a boxcar-smoothed (30 
km s^^ width), weighted average locus of emis- 
sion observed in the medium-sensitivity epochs. 
For spots < 250 mJy, the distribution is Gaussian 
and is consistent with measurement noise (Figure 
8, top panel). For spots > 250 mJy, the distribu- 
tion exhibits signiflcant non-Gaussian wings and 
tail (Figure 8, bottom panel), which is a result 
of source structure on scales of order 10 /ias. If 
we restrict the velocity range under consideration, 
the distribution of strong spots may be related 
to accretion disk thickness. The inclination and 
position angle warp of the disk combine to cre- 
ate a concavity on the near side, the bottom of 
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which is tangent to our line of sight (i.e., the disk 
is seen edge-on at this point). Maser spots in the 
velocity range 485-510 km s~^ lie in this part of 
the disk ( Herrnstein et a D I2OO5L their Figure 14). 
Previously measured accelerations for masers in 
this velocity range are 8-10 km yr^^, which 
suggests a radius of 3.7-4.1 mas and from the warp 
model we infer a projected range vertical offset of 
only 3 jias. The observed distribution of maser 
spot position offsets from the mean locus is nearly 
Gaussian with a = 5.1 /Ltas (12 /xas at full-width 
half-maximum) . Because most of the maser spot 
positions in this velocity range are measured with 
< 3/ias accuracy, we propose that this Gaussian 
distribution largely represents a true distribution 
about the local midplane of the disk (i.e., the 
thickness of warm molecular material.) See next 
section for interpretation. 

5. Notes on Source Structure and Evolu- 
tion 

The long-term importance of VLBI monitoring 
of the NGC4258 H2O maser includes contribu- 
tions to estimation of a geometric distance and 
its application to high-accuracy measurement of 
the Hubble constant. Estimation of a distance re- 
quires additional quantification of centripetal ac- 
celerations for disk material and modeling of 3- 
D disk structure and dynamics. These will be 
treated in follow-on papers. However, we can 
make several observations based on the time-series 
VLBI data alone. 

Accretion disk thickness. — In prior analy- 
ses, scatter in low-velocity maser spot positions 
(north-south) have been used to estimate an up- 
per limit on the thicknes s of th e underlying ac- 
cretion disk. iMoran et al. ( 1995 ) obtained a limit 
of ^ lO^as by focusing on strong maser features 
toward the middle of the spectrum for a VLBI 
epoch in 1994. We have combined the results 
for all 18 VLBI epochs to obtain a time-average 
sample of emission for which the width of the dis- 
tribution (Ict) is 5.1/xas. This may be the first 
directly measured thickness for an accretion disk 
around any black hole. In hydrostatic equilib- 
rium the accretion disk is expected to have a 
Gaussian distribution of density as a function of 
height with a = r Cg/v^, where r is the radius, 
is the orbital velocity, and Cs is the sound 



speed ([Frank. King, fc Raind [2002f) . Note that 
Vrj,/r = {GM/r^y/^ = uj, which is the slope of the 
impact parameter versus velocity curve of Fig- 
ure 10 and equals 2.67 x 10~^° s~^. Hence, for 
cr — 5.1 /ias or 5.5 x lO^'' cm, we obtain Cg — 1.5 
km s^^, corresponding to a tempera ture of about 
600 K (about half the limit given by iMoran et al 



(|l995h l. This temperature is consistent with the 
400-1000 K range that i s believed ne cessary for 
maser action from H2O ( Elitzur 1992f) . Further- 
more, the inferred sound speed is comparable to 
the line width of individual maser Doppler com- 
ponents seen in spectra, which is consistent with 
maser saturation. We note that in our discus- 
sion, we have assumed that the maser distribution 
samples the full thickness of the accretion disk, 
rather than a surface layer, whose depth could 
be well below that of the physical disk. This as- 
sumption is consistent w ith the disk models of 
Neufeld fc MalonevI (|l995h for radiatively efficient 



accretion. A more detailed inspection of residuals 
from models fit to maser spot distributions (e.g., 
residuals arising from an antisymmetric offset per- 
pendicular to the disk plane, caused by irradiation 
of "upper" and "lower" surfaces of the disk for 
high-velocity red- and blue-shifted masers) is nec- 
essary for further comment. We also note that the 
measured thickness places a limit on the magnetic 
pressure. The magnetic field strength must be less 
than about 100 mG, which is about the current 
limit from the non d etection of the Zeeman effect 
(|Modiaz et al.ll2005l) . 

Highly red-shifted emission. — Red-shifted emis- 
sion at velocities ^ 1460 km s^^ has been mapped 
here for the first time. Emission at 1647km s~^ 
was first reported by iModiaz et al 
a single-dish spectrum, while the 



(120051 ) in 

1566 km s~^ 



Doppler component was previously unknown. The 
detection of these very high velocity emissions en- 
ables the mid-line of the disk, along which the 
rotation curve is apparent, to be traced to about 
30% smaller radii (« 0.11 pc) than before. This 
innermost radius is now also ~ 20% smaller than 
that for the low- velocity emission. Perhaps most 
importantly, the radius of the 1647 km s~^ emis- 
sion is small enough that curvature in the disk is 
directly observable, which will enbl e more robust 
geometric modeling of the warp (cf. iMivoshi et al 
I1995I : iHerrnstein et al.Hl998[ ). 



The low- velocity emission was shown bv lHerrnstein et al 
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(j2005l ) to originate near the bottom of a concave 
depression on the near side of the disk, where our 
hne of sight is tangent to the disk. We note that in 
contrast, the 1566 kms"^ component hes near the 
topmost portion of the spine of the disk projected 
on the sky. For an "ahen" observer viewing the 
disk edge-on along a hne of sight close to our sky 
plane, the material we perceive at 1566kms~^ 
would probably be responsible for strong low- 
velocity emission. Conceivably, the intensity of 
that emission would be boosted by seed photons 
from the northern lobe of the radio jet, much as 
the low-velocity emission visible to earthbound 
observers is believe d to be boosted by the south- 
ern lobe of the jet (jHerrnstein et al.lll997r ). The 



resulting (competitive) beaming of maser energy 
along a radial path may be responsible for the low 
intensity of the 1566 km s~^ emission, for which 
the dominant gain path is orthogonal. 

Source Symmetry. — Early VLBI maps of the 
maser emission exhibited asymmetry, wherein 
the blue-shifted high-velocity emission spanned 
a truncated range of radii compared to red-shifted 
emission (|Mivoshi et al., il995). It was unclear 
whether this was a reflection of limited sensitivity 
(since the blue-shifted emission is quite weak) or 
a ramification of unmodeled disk structure. De- 
tection of emission over a broader range of radius 
(i.e., velo city) has been accomplished since then. 



Nakai et al. (1995) first observed a broader range 
via single-dish spectroscopy, detecting emission 
near —300 and — 500kms~^ one time each, in « 20 
epochs spaced oy e r w 2 00 days. In their Figure 3, 



Herrnstein et al.l (|2005l ) also show data derived 
from a VLBI detection of blue-shifted emission 
near —300 and — 500kms~^ at one epoch in 1996 
September. With new detections shown in our 
Figure 13 at « —285 km s""'^ (7 epochs over 2-3 
years) and « — 515kms~^ (2 epochs over 8 days), 
we confirm the distribution of blue-shifted high- 
velocity emission over a ~ 0.1 pc range of radius. 
In light of the very high velocity red-shifted emis- 
sion reported here (which shows that warm molec- 
ular material extends to relatively small radii), 
and assuming mirror symmetry of the maser disk, 
we hypothesize there may be still higher velocity 
blue-shifted emission that could be detected in the 
vicinity of — 600kms~^. However, any such emis- 
sion is likely to be weak, perhaps on the order of 
1 mjy, based on line ratios among known Doppler 



components. 

Envelope of low-velocity emission. — Our long 
term monitoring of velocities between the high 
and low-velocity emission complexes has not re- 
sulted in detected intermediate-velocity emis- 
sion. Nondetection here and in episodic single- 
dish spectroscopy o f comparable sensitivity (e.g., 
Modjaz et al. 2005h may exclude the possibility of 



long gain paths at disk azimuth angles more than 
« 0.1 rad from the emission close to the systemic 
velocity. This contrasts with what has been found 
in th e NGC 1068 H2O maser ( GreenhiU fc Gwinn 
Il997h . wherein the locus of "low-velocity" maser 
emission extends in a ring over a full quadrant 
on the receding side of the disk. However, in 
the angular extent of high brightness temperature 
continuum emission that can provide seed pho- 
tons for amplification differs. In NGC 4258, they 



arise from a compact jet core (jHerrnstein et al 



19971 . 120051 ). while in NGC 1068, they are believed 



to arise from an extended region that more or 
less fi lls the volume inside the innermost maser 
orbit (jGallimore et al. I l200lh -though the absence 
of seeded emission in the adjacent (approaching) 
quadrant is difficult to explain. 

Though centripetal acceleration is responsi- 
ble for the observed secular drift of individual 
Doppler component velocities, the overall range 
of emission does not change discernably in angle 
or velocity. In our time monitoring, we detected 
low- velocity emission from w 380 to 580kms~^. 
This range is comparable to the ranges ob- 
served in high-sensitivity spectr oscopic studies 



over t h e period 1992-2005 (e.s.. I GreenhiU et al 
1995bl: iNakai et all Il995l : lYamauchi et al.l l2005t 



Modiaz et al.l iooiE . Any secular change of the 



upper and lower limits of the envelope appears to 
be much smaller than the lOkms^^yr^^ drift of 
individual Doppler components. 

No maser spots seen away from the disk. — 
In contrast to maser emission f rom the Circi- 
nus galaxy d GreenhiU et ai~1 l2003l) and NGC 1068 
( Gallimore et al.l 19961 2001 ) , we found no Doppler 
components more than a beamwidth (w 1 mas) 
from the locus of disk emission. Our detection 
limit was 13 mJy (6cr) for the highest sensitivity 
observation, BM112H. In Circinus, the maser pop- 
ulation not associated with the disk is believed to 
trace a bipolar, wide-angle outfiow that contains 
clumps moving away from the central engine, at 



11 



radii as small as « 0.1 pc. In NGC 1068, the 
non-disk emission lies ^ 1 pc downstream in the 
narrow- line region, along the jet axis. The absence 
of emission associated with wind or jet emission 
in NGC 4258 (see §2.4) may be a consequence of a 
lower luminosity central engine and more tightly 
bound disk material (i.e., smaller disk mass com- 
pared to the central engine mass). 
Absorption signatures. 



Watson k Wallinl(|l994^ 



hypothesized that thermalized water molecules 
in disks could absorb maser emission arising at 
smaller radii; the observational signature would 
be troughs in spectra near systemic velocities. 
A number of spectra for NGC 4258 that are in 
the literature do exhibit dips in the vicinity of 
the systemic velocity, including time averages of 
many spectra collected over pe riods of months 



or years (see iBragg et al.l l2000l and references 



therein, as well as our Figures 1,5). The cumula- 
tive record of VLBI observations does not demon- 
strate notable structure in the vicinity of the sys- 
temic velocity, which is consistent with expecta- 
tions. However, evidence thus far has been anec- 
dotal (e.g., compare spectra for epochs 2000.08 
and 2000.61). Systematic analysis of spectra ob- 
tained throughout 1985 and 1986 demonstrated 
that troughs in spectra exhibited secular drift 
similar to that of neighbori ng emission compo- 
nents (jGreenhill et al.lll995bh . A physical model 
of the NGC 4258 warped accretion disk presents 
additional hurdles to the hypothesis. In the case 
of NGC 4258, the vertical deviation of the disk 
from a plane disp laces material at large radii from 
the line of sight ( Herrnstein et al. 20051 ). and the 
velocity at which the longest coherence lengths 
occur (overlapping the rest frame transitions of 
inverted and thermalized material) is displaced 
^ 20kms^^ redward from the systemic velocifr 
Moreover, the model of .Ncufcld_&LMaloney (199 
predicts that material at sufficiently large radii is 
atomic, and thus cannot form an absorbing layer. 

6. Conclusion 

We have reported a four year time-series VLBI 
study of the II2O maser in NGC 4258, with 
quadruple the number of epochs than earlier work. 
Overall source structure during the four years of 
our study resembled what has previously been fit- 
ted to an inclined, warped, Keplerian disk model. 



Frequent sampling over years enabled compila- 
tion of a more complete census of time variable 
Doppler components and broad observing band- 
width enabled detection and mapping of material 
with larger orbital velocities than before. The 
persistence and symmetry of source structure is 
now more clear, and the accuracy of position mea- 
surement has been improved. Detailed analysis of 
the distribution of maser spots perpendicular to 
the local disk plane, interpreted in the context of 
hydrostatic support, suggests a disk thickness {a 
— 5.5 X 10^* cm) and gas temperature of about 
600 K. The new catalog of maser spot positions 
and line-of-sight velocities will enable construction 
of a more robust (3-D) dynamical disk model and 
more accurate measurement of maser centripetal 
acceleration and proper motion. Together, these 
will enable estimation of a geometric distance 
wit h smaller randorn and systematic uncertainty 
(cf. iHerrnstein et al. 19991 ) . Follow-up papers will 
treat the measurement of acceleration, modeling, 
and comparison of a final "maser distance" to a 
new four-color Cepheid distance that enables re- 
calibration of distance indicators and ultimately 
Ho. 
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Table 1 
The Observations 



Experiment Code 




Date 


Antennas* 


Synthesized Beam^ 
(mas X mas, deg) 


Sensitivity*^ 
(mJy) 


O 

L 


overa^ 
R 


B 


Comments 


BM056C 


1997 


March 06 (1997.18) 


VLBA, VLA, EFLS 


0.65x0.53, 


12.5 


3.0 


X 


X 


X 




BM081A 


1997 


October 01 (1997.75) 


VLBA, VLA, EFLS 


0.59x0.35, 


-7.2 


2.7 


X 


X 


X 




BM081B 


1998 


January 27 (1998.07) 


VLBA, VLA, EFLS 


0.60x0.42, 


5.1 


3.2 


X 


X 


X 




BM112A 


1998 


September 05 (1998.68) 


VLBA, VLA, EFLS 


0.67x0.49, 


8.5 


3.0 


X 


X 


X 




BM112B 


1998 


October 18 (1998.80) 


VLBA 


0.52x0.36^ 


-18.0 


4.2 


X 


X 






BM112C 


1998 


November 16 (1998.88) 


VLBA 


0.53x0.35, 


-14.9 


4.0 


X 




X 




BM112D 


1998 


December 24 


VLBA 








X 


X 




miscorrelated 


BM112E 


1999 


January 28 (1999.08) 


VLBA 


0.51x0.34, 


-6.0 


4.4 


X 




X 




BM112F 


1999 


March 19 (1999.21) 


VLBA 


0.48x0.34, 


-10.0 


4.1 


X 


X 






BM112G 


1999 


May 18 (1999.38) 


VLBA 


0.50x0.35, 


-10.8 


4.8 


X 




X 




BM112H 


1999 


May 26 (1999.40) 


VLBA, VLA, EFLS 


0.43x0.26, 


-31.3 


2.3 


X 


X 


X 




BM112I 


1999 


July 15 


VLBA 








X 


X 




corrupted 


BM112J 


1999 


September 15 (1999.71) 


VLBA 


0.51x0.32, 


-5.4 


5.8 


X 




X 




BM112K 


1999 


October 29 (1999.83) 


VLBA 


0.50x0.34, 


-19.2 


4.3 


X 


X 






BM112L 


2000 


January 07 (2000.02) 


VLBA 


0.50x0.34, 


-12.0 


3.6 


X 




X 




BM112M 


2000 


January 30 (2000.08) 


VLBA 


1.25x0.35, 


38.8 


4.3 


X 


X 






BM112N 


2000 


March 04 (2000.17) 


VLBA 


0.48x0.33, 


-10.9 


3.9 


X 




X 




BM1120 


2000 


April 12 (2000.28) 


VLBA 


0.63x0.37, 


-24.4 


5.0 


X 


X 






BM112P 


2000 


May 04 (2000.34) 


VLBA 


0.50x0.34, 


-11.4 


4.6 


X 




X 




BG107 


2000 


August 12 (2000.61) 


VLBA, VLA, EFLS 


0.46x0.36, 


-14.3 


4.7 


X 


X 


X 





^VLBA: Very Long Baseline Array; VLA: 27x25-m Very Large Array in Socorro, NM; EFLS: Max-Planck- Inst it ut fur Radioastronomie 100-m 

antenna in EfTelesberg, Germany 

^Average restoring beam full-width half-power and position angle, measured east of north. The unusually large synthesized beam in BM112M 
is due to the loss of two antennas, HN and NL. See data tables for more detail. 

*^ Average rms noise (scaled to a lkms~^ channel spacing) in emission-free portions of synthesis images. See Table 5 for more detail. 
'^Velocity range observed. S: low-velocity emission; R: red-shifted emission; B: blue-shifted emission 



Table 2 
Instrument Setup 



Experiment Code Bandwidth Emission Type 


Bandpass 


c 

Band center vlsr 


• 

Tuning of Bandpass 


(MHz) 




/I — 1 \ 
(km s ) 




DiviuoDk-') iJiviuoi/i.j Divii izi\,xi o low-veiocicy 


]^ 


1 s so 


LCP 




2 


ooo.yo 


LCP 




q 


511 24 


LCP 






426 30 


LCP 


red-shifted 


5 


1447.41 


RCP 




Q 


1 S^Q S7 


RCP 




7 


1339 76 


RCP 






1251 72 


RCP 


Diue-snirteci 


9 


— 33fi flQ 


RCP 




10 


—424.00 


RCP 




11 


—443.74 


RCP 




12 


— 531 65 


RCP 


BG107 8 low-velocity 


1 


665.00 


LCP 




2 


580.00 


LCP 




3 


495.00 


LCP 






/tin nn 

^lU.UU 


LCP 


red-shifted 


5 


1485.00 


RCP 




6 


1400.00 


RCP 




7 


1315.00 


RCP 




8 


1230.00 


RCP 


blue-shifted 


9 


-255.00 


RCP 




10 


-340.00 


RCP 




11 


-425.00 


RCP 




12 


-510.00 


RCP 


BM112B,F,K,M,0 16 low-velocity 


1 


836.00 


LCP 




2 


653.00 


LCP 




3 


470.00 


RCP 




4 


470.00 


LCP 


red-shifted 


5 


1568.00 


LCP 




6 


1385.00 


LCP 




7 


1202.00 


LCP 




8 


1019.00 


LCP 


BM112C,E,G,J,L,N,P 16 low velocity 


1 


500.00 


RCP 




2 


500.00 


LCP 




3 


317.00 


LCP 




4 


134.00 


LCP 


blue-shifted 


5 


-49.00 


LCP 




6 


-232.00 


LCP 




7 


-415.00 


LCP 




8 


-598.00 


LCP 


^See Table 1 for dates of observation. 








''The 8 MHz and 16 MHz bands span 108 km s~^ and 216 km s~\ 


respectively. 






^Velocities are classical radio velocities in the Local Standard of Rest (LSR) reference frame. 
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Table 3 

Station Positions and Velocities'* 



Station Station Position Error 



X Y Z X Y Z 





(m) 


(m) 




(m) 




(cm) 


( 


cm) 


(■ 


cm) 


vlba-br 


-2112064.9661 


-3705356, 


.5115 


4726813, 


.7932 


0.0476 


0, 


,0885 


0, 


1069 


EFLS 


4033947.4792 


486990 


.5244 


4900430, 


,8031 


0.2577 


0, 


,0669 


0, 


,3094 


VLBA-FD 


-1324009.1202 


-5332181, 


.9713 


3231962 


,4738 


0.0327 


0, 


1031 


0, 


,0731 


VLBA-HN 


1446375.1259 


-4447939, 


,6532 


4322306, 


,1181 


0.0395 


0, 


1275 


0, 


1237 


VLBA-KP 


-1995678.6199 


-5037317 


,7147 


3357328, 


,1296 


0.0453 


0, 


1104 


0, 


,0832 


VLBA-LA 


-1449752.3520 


-4975298, 


,5907 


3709123, 


,9270 


0.0290 


0, 


,0830 


0, 


,0689 


VLBA-MK 


-5464074.9591 


-2495249, 


.1189 


2148296, 


.8444 


0.1700 


0, 


,1071 


0, 


,1037 


VLBA-NL 


-130872.2452 


-4762317, 


.1215 


4226851, 


.0426 


0.0217 


0, 


,1067 


0, 


,0953 


VLBA-OV 


-2409150.1035 


-4478573, 


.2283 


3838617, 


.3949 


0.0546 


0, 


,1027 


0, 


,0919 


VLBA-PT 


-1640953.7026 


-5014816, 


.0276 


3575411, 


.8817 


0.0321 


0, 


,0851 


0, 


,0692 


VLBA-SC 


2607848.5236 


-5488069, 


.6773 


1932739, 


.5281 


0.1088 


0, 


,2388 


0, 


,1104 


VLA 


-1601185.2077 


-5041977, 


.1729 


3554875, 


.7033 


fs 1.0 




a 1.0 




a 1.0 



Station Station Velocity Error 



X Y Z X Y Z 

(cm yr~^) (cm yr""*^) (cm yr~^) (cm yr""*^) (cm yr""*^) (cm yr~^) 



VLBA-BR 


-1.2800 


0.0370 


-1.1200 


0.0160 


0.0280 


0.0304 


EFLS .... 


-1.5250 


1.7100 


0.8950 


0.0492 


0.0184 


0.0626 


VLBA-FD 


-1.1600 


-0.2380 


-0.7040 


0.0102 


0.0234 


0.0194 


VLBA-HN 


-1.4430 


-0.1430 


0.0850 


0.0142 


0.0436 


0.0387 


VLBA-KP 


-1.2210 


0.0150 


-0.9740 


0.0149 


0.0309 


0.0255 


VLBA-LA 


-1.2850 


-0.0170 


-0.8300 


0.0083 


0.0174 


0.0165 


VLBA-MK 


-1.3600 


6.0080 


2.9420 


0.0547 


0.0428 


0.0427 


VLBA-NL 


-1.4180 


0.0090 


-0.4710 


0.0095 


0.0289 


0.0245 


VLBA-OV 


-1.6270 


0.5910 


-0.8200 


0.0169 


0.0305 


0.0268 


VLBA-PT 


-1.3610 


-0.2480 


-0.8590 


0.0073 


0.0147 


0.0151 


VLBA-SC 


0.9840 


0.4390 


1.0200 


0.0331 


0.0784 


0.0445 


VLA 


-1.3140 


-0.1100 


-0.8560 


Ri 0.1 


Ri 0.1 


« 0.1 



^The VLBA and EFLS station positions and velocities are from the United States Naval Observatory's 
(USNO's) VLBI n98 10 solution, which has a reference date of January 1, 1997. USNO's global VLBI solutions 
determine the International Terrestrial Reference Frame (ITRF), the International Celestial Reference Frame 
(ICRF), and Earth Orientation Parameters (EOPs). These solutions, therefore, yield a consistent set of station 
positions and velocities, source positions, and Earth Orientation Parameters. The VLA station position and 
velocity is from an NRAO 12/00 geodesy solution (Craig Walker, private communication). It has been rotated 
to the ITRF reference frame and also has a reference date of 1997 January 1. 
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Table 4 

Largest Sources of Relative Position Error in Synthesis Maps'* 



Type of Error 


Magnitude over 75 MHz'' 
(fias) 


Magnitude over 4 MHz'' 
(^as) 


Comments 


Random^ 


1 
25 


1 
25 


1 Jy peak (strong) 
50 mjy peak (weak) 


Systematic (full frequency range): 








Positional Uncertainty of Reference Feature^ 
Atmospheric Delay Error* 


10 
6 


< 1 

< 1 


Uncertainty in absolute position of reference feature: 
a = 0f0003, S = o:'003 


Frequency Standard Uncertainty® 


2 


< 1 




Positional Uncertainty of Calibrator*^ 
Baseline Coordinate Uncertainty^ 


1 
1 


< 1 

< 1 


Uncertainty in absolute position of calibrators: 
a = 0f000042, 6 = 0"00047 (4C39.25), 
a = 0f000114, S = 0"00027 (1150+81), 
a = 0f000020, S = 0"00026 (1308+33) 



Systematic (relative frequencies): 



Clock Uncertainty^ ... 2 



^The errors discussed in this table arc analytic and affect all epochs. A few epochs are subject to an additional, non-analytic source of systematic error, which is 
discussed in §3.2 of the text and in Table 7. 

"^75 MHz is the characteristic frequency separation between low- velocity and red-shifted or low- velocity and blue-shifted masers. Since everything is referenced to 
a low-velocity maser, this is the maximum frequency separation that need be considered. 

'^The full 75 MHz frequency span is applicable to single epoch position offsets only. When comparing one epoch's position offsets to another's (as with proper 
motions), a 4 MHz frequency span is more realistic. This is because the reference channels for the epochs in Table 1 span only 4 MHz, meaning that systematic 
errors largely cancel for red- and blue-shifted masers. 

■^The formal uncertainty in the fitted position of an unblended maser spot in a single channel is given by AG^^t ~ 0.bO-\-,{AS/ S), where is the synthesized beam 
size, AS is the image rms, and S is the peak intensity JReid et al.||198^ . The table values assume — 500 fias and AS — 5 mJy. Random errors are not frequency 
dependent and do not cancel when comparing position offsets of a given maser feature between epochs. 

^The relative position error in synthesis maps due to the uncertainty in the absolute position of the reference feature is given by (Ai//f) A^ref , where Ai^ — u — t'rcf , 
i^ref is the reference frequency, and A^rcf is the uncertainty in the absolute position of the reference feature (jThompson. Moran, Swensonl|200in . 

^The relative position error in synthesis maps due to the atmospheric delay error remaining after delay calibration corrections have been applied to the maser 
sourc e is given by (c/B) [Au/i') Sta-, where c is the speed of light, B is the baseline length, and 5ta is the atmospheric delay error l lThompson. Moran, Sz SwensonI 
120011 ). Other parameters are defined above. 5TA{'nsec) ~ sec z tan z Az, where z is the zenith angle and Az is the difference in zenith angle between the delay 
calibrator used to make the corrections and the maser source. Taking average maser-calibrator offsets and zenith angles of ^ 15° and ~ 35°, respectively, we get a 



6 fias relative position error for features 75 MHz from the reference feature on an 8000 km baseline, appropriate for EFLS— VLA, the most sensitive baseline. 



^In general, errors in station clocks consist of two components: a fixed offset and a time dependent term arising from frequency standard errors. Hydr ogen maser 
frequency standards typically have a fractional uncertainty of 10~^^, which leads to a residual fringe rate of 1 mHz. iThompson, Moran. fc SwensonI ( |200lD show that 
the resulting relative position error in synthesis maps is given by {c/B) (A.iy/v) {6{Si^) /lyo) t, where 5(5f)/t'o is the fractional fringe rate uncertainty with ivq being 
the LO frequency, and t is the time between calibrator scans. Other parameters are defined above. In computing table values, we assumed a 20-minute calibrator 
separation and an 8000 km baseline. 

^An uncertainty in absolute calibrator positi on leads to a residual fringe rate (in mHz) of 0.13 {B/b X 10'^ km) mas), where B is the base li ne length 

and AOcai is the position error of the calibrator ^Thompson, Moran, &: Swensonll 20011) . When applied to the line source, iThompson, Moran, Sz SwensonI < i26'oil) show 
that this leads to relative position errors in synthesis maps with the same functional dependence as above, but with t being the total integration time rather than 
the time between calibrator scans. It is difficult to assess the combined error of the multiple calibrators and many baselines in synthesis maps, but it is likely that 
the dominant contribution in most high sensitivity epochs comes from an approximately 2.5 hour 4C39.25 time span on the EFLS— VLA baseline, which is the most 
sensitive. We note that calibrator coordinate errors are given with respect to the International Celestial Reference Frame (ICRF) and are consistent with the Earth 
Orientation Parameters (EOPs) used during correlation and subsequent data analysis. 

'The relative positio n error in synthesis maps due to a ba seline error is given by (Ai'/i-') {AB / B), where is defined as above, AB is the baseline error, and B 
is the baseline length ^Thompson, Moran, & Swenso iil l2001D . The EFLS-VLA baseline is representative, since, being the most sensitive, it is expected to have the 
largest effect on synthesis maps. As mentioned in the text, the position of the VLA is uncertain to 3 cm, while all other station positions are uncertain to 1 mm. 
Baseline coordinates are consistent with the ICRF and the EOPs (see note above). 

■•Clock offset error s do not scale with the f ull frequency range, but rather with the relative frequency of remote and reference channels, each with respect to their 
own LO freouencies (]Thompson, Moran, &z Sw enson 20 ^). Hence, Ai/j-ei — (i-* — fo) — (j-'re/ — ^T-e/o ) ■ The relative position error in synthesis maps due to clock errors 
is then {Avrei / ^){cSto / B) , where Stq is the post calibration delay residual. The accuracy of delay measurements is expected to be of or der (l/2Ai/(,u;) (l/SNR) , 
where Afb^^ is the total bandwidth and SNR is the signal-to-noise ratio of calibrator measurements ^Thompson. Moran, & Swenson||200lD . Assuming SNR — 150, 
we obtain a residual delay of 0.4 nsec and a relative position error in synthesis maps of 2 fias for relative frequency separations of 4 MHz. 



Table 5 
Channel Data 



channel ^^lsr'^ 








Ax" 




Ay" 


c -d d-d ddd c 
<7y maj CTmaj mm (Tmin pa CTpa ^^band 


(1-512) (kms^^) 


(km s"^) 


(Jy) 


(mJy) 


(mas) 


(mas) 


(mas) 


(mas) (mas) (mas) (mas) (mas) (°) (°) (kms^^) 



177 


1464, 


.157 


1467, 


702 


0, 


.038 


3 


.690 


4, 


.458 


0, 


.024 


0. 


,784 


0, 


.032 


0. 


82 


0, 


,08 


0, 


.57 


0, 


.05 


167 


10 


1447. 


,41 


178 


1463, 


.946 


1467. 


,490 


0, 


.054 


3, 


.850 


4, 


.486 


0, 


.018 


0. 


,783 


0, 


.023 


0. 


,71 


0. 


,05 


0, 


.53 


0, 


.04 


7 


10 


1447. 


,41 


179 


1463, 


.735 


1467. 


,279 


0, 


.057 


3, 


.770 


4, 


.499 


0, 


.017 


0. 


,750 


0, 


.022 


0. 


,71 


0. 


,05 


0, 


.54 


0, 


.04 


19 


9 


1447. 


,41 


180 


1463, 


.525 


1467. 


,067 


0, 


.044 


3, 


.680 


4, 


.508 


0, 


.021 


0. 


,753 


0, 


.028 


0. 


68 


0. 


,06 


0, 


.60 


0, 


.05 


23 


27 


1447. 


,41 


223 


1454, 


.466 


1457. 


,964 


0, 


.027 


3, 


.570 


4, 


.475 


0, 


.033 


0. 


,778 


0, 


.043 


0. 


,71 


0. 


,09 


0, 


.42 


0, 


.06 


12 


10 


1447. 


,41 


224 


1454, 


.255 


1457. 


,752 


0, 


.039 


3, 


.550 


4, 


.503 


0, 


.023 


0. 


,791 


0, 


.030 


0. 


65 


0. 


,06 


0, 


.49 


0, 


.04 


14 


13 


1447. 


,41 


225 


1454, 


.045 


1457. 


,541 


0, 


.057 


3, 


.560 


4, 


.531 


0, 


.016 


0. 


,776 


0, 


.020 


0. 


65 


0. 


,04 


0, 


.55 


0, 


.03 


14 


15 


1447. 


,41 


226 


1453, 


.834 


1457. 


,329 


0, 


.079 


3 


.620 


4, 


.523 


0, 


.012 


0. 


,782 


0, 


.015 


0. 


68 


0. 


,03 


0, 


.55 


0, 


.03 


24 


8 


1447. 


,41 


227 


1453, 


.623 


1457. 


,117 


0, 


.092 


3, 


.730 


4, 


.513 


0, 


.010 


0. 


,800 


0, 


.013 


0. 


,71 


0. 


,03 


0, 


.55 


0, 


.02 


15 


6 


1447. 


,41 


228 


1453, 


.413 


1456. 


,906 


0, 


.102 


3 


.720 


4, 


.512 


0, 


.009 


0. 


,814 


0, 


.012 


0. 


69 


0. 


,03 


0, 


.55 


0, 


.02 


9 


6 


1447. 


,41 


229 


1453, 


.202 


1456. 


,694 


0, 


.106 


3, 


.620 


4, 


.515 


0, 


.009 


0. 


,802 


0, 


.011 


0. 


68 


0. 


,02 


0, 


.53 


0, 


.02 


8 


5 


1447. 


,41 


230 


1452, 


.991 


1456. 


,482 


0, 


.112 


3, 


.610 


4, 


.517 


0, 


.008 


0. 


,803 


0, 


.011 


0. 


65 


0. 


,02 


0, 


.52 


0, 


.02 


8 


6 


1447. 


,41 


231 


1452, 


.781 


1456. 


,271 


0, 


.114 


3 


.560 


4, 


.515 


0, 


.008 


0. 


,818 


0, 


.010 


0. 


65 


0. 


,02 


0, 


.54 


0, 


.02 


6 


6 


1447. 


,41 


232 


1452, 


.570 


1456. 


,059 


0, 


.115 


3, 


.680 


4, 


.538 


0, 


.008 


0. 


,823 


0, 


.011 


0. 


66 


0. 


,02 


0, 


.55 


0, 


.02 


17 


7 


1447. 


,41 


233 


1452, 


.359 


1455. 


,847 


0, 


.116 


3, 


.680 


4, 


.546 


0, 


.008 


0. 


,821 


0, 


.010 


0. 


69 


0, 


,02 


0, 


.52 


0, 


.02 


9 


4 


1447. 


,41 



Note. — Tabic 5 provides the Channel data for all the observations and is published in its entirety in the cleetronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 

^Velocities according to the classical and rclativistic definitions of the Dopplcr shift in the LSR frame. Sec §2.6 for a detailed discussion of conversion 
from classical to r elativistic LSR velocity. Classical Dopplcr tracking was implemented using the AIPS task CVEL, which is accurate to < 0.004 km s~^ 
JBragg et al.|[2000i) . 

^Fitted peak intensity (Jy beam"^) is reported in this table. Note that Flux density (Jy) — Fitted peak intensity (Jy bcam"^) ^{Ol^^cis/^bea-m)^ whcrc 
^meas mcasurcd maser spot area and ^beam area of the restoring beam. Since maser spots arc assumed to be unresolved, these values can be 

quoted in Jy. 

'^East-west and north-south position offsets and uncertainties, measured with respect to the emission at 510.0 km s^^ (relativistic). Position uncertainties 
reflect fitted random error only. The true random error is estimated to be the maximum of the fitted random error and the random error computed according 
to Table 4, footnote d. 

^Major axes, minor axes, and position angles of maser spots and uncertainties. 

^Band identifier refers to the nonrelativistic velocity in channel 256.5. 



Table 6 
Continuum Detections 



Epoch'^ 


Freq.'' 




rms^ 






Ay" 


ay 






(mJy) 


(mJy) 


(mas) 


(mas) 


(mas) 


(mas) 


BM056C 


blue-shifted 


1.19 


0.18 


-0.240 


0.051 


0.952 


0.055 


BM081B 


red-shifted 


1.66 


0.19 


0.108 


0.022 


0.800 


0.032 




blue-shifted 


1.28 


0.19 


-0.194 


0.028 


1.018 


0.035 


BM112B 


red-shifted 


1.53 


0.19 


0.007 


0.037 


1.360 


0.099 


BM112C 


blue-shifted 


1.67 


0.18 


-0.190 


0.028 


1.099 


0.070 


BM112E 


blue-shifted 


2.40 


0.20 


-0.028 


0.030 


0.952 


0.048 


BM112F 


red-shifted 


1.90 


0.19 


-0.039 


0.054 


0.897 


0.093 


BM112G 


blue-shifted 


0.72 


0.21 


0.200 


0.047 


1.008 


0.080 


BM112J 


blue-shifted 


1.70 


0.26 


-0.147 


0.075 


1.093 


0.147 


BM112K 


red-shifted 


1.72 


0.10 


-0.192 


0.134 


1.534 


0.213 


BM112L 


blue-shifted 


1.75 


0.16 


-0.009 


0.029 


1.252 


0.061 


BM112M 


red-shifted 


2.68 


0.20 


-0.346 


0.363 


0.990 


0.157 


BM112N 


blue-shifted 


2.39 


0.18 


-0.049 


0.013 


0.871 


0.030 


BM1120 


red-shifted 


4.01 


0.24 


0.022 


0.021 


0.947 


0.040 


BM112P 


blue-shifted 


2.52 


0.21 


-0.029 


0.061 


1.018 


0.077 


BG107 


red-shifted 


3.95 


0.22 


-0.133 


0.015 


1.050 


0.018 




blue-shifted 


3.66 


0.20 


-0.152 


0.010 


1.022 


0.013 



^Wc report epochs where continuum emission is observed at > 3cr 

'^For BM056C and BM081B we averaged the inner 490 channels of the three inde- 
pendent red-shifted and three independent blue-shifted bands (1470 channels). For 
all other epochs we averaged the inner 470 channels of the four red-shifted and four 
blue-shifted bands (1960 channels). 

'^Wc report the flux density, since coniiniium detections are partially resolved. Flux 
density (mJy) — Fitted peak intensity (mJy beam""'") X (^^eas /^beam )> where O'^^^^ 
is the measured area of the continuum emission and O^^am area of the restoring 

beam. Strictly speaking, this is not a total flux density, but rather an angle integrated 
flux density, since we measure it with an interferometer that has a limited range of 
baseline lengths. 

^la noise level in synthesis map. 

^East-west and north-south position offsets and uncertainties measured with re- 
spect to the reference maser feature adopted for the epoch of observation. Position 
uncertainties reflect random error only. 
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Table 7 
Derr'ed Positiox Shifts 



Experiment Code 


L. 




In- Mas, 




1 


)<-il-slnH< 


•<1 Masr 




X 

(/xas) 


(^as) 


y 

(/xas) 


(/xas) 


X 

(/A as) 


(/xas) 


y 

(/A as) 


(Ty 

(/A as) 


BM056C 


8 


3 


1 


4 


5 


5 


-1 


6 


BM081A 


5 


1 


-3 


2 


-8 


4 


-10 


6 


BM081B 


-4 


2 


-11 


2 


-1 


3 


-3 


5 


BM112A 


-8 


2 


3 


4 


33 


5 


6 


8 


BM112B 


5 


3 





5 


2 


5 





7 


BM112C 


3 


2 


-2 


3 










BM112E 





3 


1 


4 










BM112F 


4 


2 


-2 


3 


-3 


4 


-6 


6 


BM112G 


-3 


3 


8 


4 










BM112H 


-6 


1 





2 


-24 


4 





6 


BM112J 


10 


3 


-6 


4 










BM112K 


1 


2 


-1 


3 


-5 


6 


3 


8 


BM112L 


3 


2 


-3 


3 










BM112M 


-7 


3 


-5 


10 


-5 


4 





14 


BM112N 


3 


3 


6 


4 










BM1120 


-7 


4 


9 


6 


-12 


4 


15 


7 


BM112P 


-2 


3 


-10 


4 










BG107 


-20 


4 





5 


-17 


7 


37 


9 



^Position shifts are determined by minimizing the sum of the deviations in each epoch from 
the weighted average of the masers in the medium sensitivity epochs. Weighted averages 
are computed and deviations determined for each frequency channel. See §3.2 for details. 
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Fig. 1. — Spectra for fow-velocity masers. The 

vertical scale on the page is the same for each 
plot. The spectra show the total imaged power ob- 
tained from the fitted peak flux densities of maser 
spots. Channels with no detectable emission are 
not shown. Velocities are defined with respect to 
the Local Standard of Rest, assuming the nonrel- 
ativistic radio definition of Dopplcr shift (bottom 
axis) and relativistic radio definition of Doppler 
shift (top axis). 
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Fig. 2. — Spectra for rcd-shiftcd masers, dis- 
played as in Figure 1. The Doppler compo- 
nents at 1566kms-i (epochs 1999.21, 1999.83, 
and 2000.08) and 1647 km s'^ (epochs 1999.83, 
2000.08, and 2000.28) are shown separately in Fig- 
ure 3. 
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Relativistic LSR Velocity (km s ^) 




Radio LSR Velocity (km s^^) 

Fig. 3. — Spectra for red-shifted masers at 
1566 kms~^ and 1647kms~^, displayed as in Fig- 
ure 1. The dotted portions of the spectra deUneate 
fits that fall below our 6a threshold (see text for 
details). 
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Relotiviatic LSR Velocity (km a '') 




Fig. 4. — Spectra for blue-shifted masers, dis- 
played as in Figure 1. The dotted portions of the 
spectra delineate fits that fall below our 6a detec- 
tion threshold (see text for details). 
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Relativistic LSR Velocity (km s ^) 
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Fig. 5. — Unweighted average maser spectra for 
all observing epochs (1997-2000), displayeci as in 
Figure 1. (top) Low- velocity emission, (middle) 
Red-shifted high- velocity emission, (bottom ) Blue- 
shifted high-velocity emission. LSR velocities are 
defined assuming the nonrelativistic radio defini- 
tion of Doppler shift (bottom axis) or relatlvistic 
radio definition of Doppler shift (top axis) . 
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Fig. 6. — Maps of maser emission for all epochs 
superposed, (top) The full velocity range of the 
emission, (upper middle) Low- velocity emission. 
(lower middle) Red-shifted high- velocity emis- 
sion, (bottom) Blue-shifted high- velocity emission. 
Symbol diameter is proportional to the logarithm 
of peak flux density. Each panel is scaled sepa- 
rately. A key on the right hand side of the three 
lower panels gives the symbol sizes of 2.5 Jy, 250 
mjy, and 25 mJy masers, as well as the x-y sizes of 
their associated random errors. The color bar at 
bottom codes nonrelativistic LSR velocity (lower 
label) or relativistic LSR velocity (upper label). 
Note that the thinness of the maser disk in the 
low- velocity panel is obscured by the large number 
of masers plotted. See Figures 7-9 for additional 
discussion. 
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Fig. 7. — Map of low-velocity maser emission. 
All masers have the same symbol size. The red 
curve is the boxcar smoothed (30 km width) 
weighted average of low-velocity masers, com- 
puted from medium-sensitivity epochs (see §3.2). 
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Fig. 8. — Histograms of separations between in- 
dividual low- velocity maser spots and the average 
emission locus (Figure 7, red curve). Each off- 
set is divided by the corresponding random mea- 
surement uncertainty. All low-velocity masers are 
plotted, (top) Low- velocity masers < 250 mJy, 
for which corresponding measurement uncertain- 
ties are > 5/xas, plotted. The distribution for 
weak masers is Gaussian (see superposed curve), 
as would be expected from a noise dominated dis- 
tribution, (bottom) Low- velocity masers > 250 
mJy, for which corresponding measurement uncer- 
tainties are < 5|Uas, plotted. The distribution for 
strong masers exhibits nonGaussian wings and tail 
(see superposed curve) that may be attributed to 
source structure or systematic error arising from 
the representation of the disk plane position by 
the average of maser spot positions (Figure 7, red 
curve), at the few fj.as level. 
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Fig. 9. — Histogram as in Figure 8 for masers 
> 250 mJy but for the velocity range 485 — 510 
km s"^ (East Offset -0.1 to 0.0 mas). Tliis is 
the veloci ty range over which t he model accre- 
tion disk ( Herrnstein et al. 20051 ) is viewed edge- 
on. The distribution of strong masers is nearly 
Gaussian (see superposed curve) with full- width at 
half-maximum (FWHM) of ~ 12/ias (cr ~ S.l^as), 
which we interpret to be the underlying accretion 
disk thickness. 
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Fig. 10. — Position-velocity diagram for all epochs 
superposed. Positive impact parameters are to 
the east. The linear gradient in low-velocity emis- 
sion is characteristic of emission from material in 
highly inclined circular orbits over a narrow range 
of radii. The fitted straight line marks the locus 
for emission at a single radius. (A steeper line 
indicates a smaller radius). High- velocity masers 
trace declining rotation curves, which are fitted 
here in position-velocity space and shown for Ke- 
plerian orbits. Good agreement with the model 
V — Vsys oc where Vsys is systemic veloc- 

ity and b is impact parameter, suggests that the 
high- velocity emission arises close to a single diam- 
eter through the underlying disk an d that warping 



is of second-order importance (e.g.. iMivoshi et al 



19951 ). We note that the Doppler components at 
1566 km s^^ and 1647 km s^^ appear to arise in- 
side the mean radius of the low-velocity emission. 
Nothing has previously been mapped at these 
small radii. LSR velocities are defined assuming 
the nonrclativistic radio definition of Doppler shift 
(left-hand axis) or relativistic radio definition of 
Doppler shift (right-hand axis). 



32 



600 - 



550 



500 



o 
o 
1) 

> 
a: 



X) 
D 

d; 



450 



400 




- 600.6 



550.5 



a 



500.4 I 



Ln 

< 

o 
o 



450.3 



400.3 



-0.5 



Impact Parameter 



Fig. 11. — Position-velocity traces for low- velocity 
emission at each epoch individually. The impact 
parameter scale is shown in the bottom left corner 
of the figure and the separation between traces is 
proportional to the time between epochs (denoted 
by numeric labels). Horizontal error bars denote 
uncertainty in impact parameter and reflect ran- 
dom error only. LSR velocities are defined assum- 
ing the nonrclativistic radio definition of Doppler 
shift (left-hand axis) or relativistic radio definition 
of Doppler shift (right-hand axis). 
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Fig. 12. — Position-velocity traces for high- 
velocity red-shifted emission, plotted as in Fig- 
ure 11. There arc fewer traces here because red- 
shifted emission was only observed in alternate 
medium-sensitivity epochs (Table 1). Keplerian 
curves (sec Figure 10) were drawn through the 
data of each epoch for ease of display. 
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Fig. 13. — Position-velocity traces for high- 
velocity blue-shifted emission, plotted as in Fig- 
ure 11. There arc fewer traces here because blue- 
shifted emission was only observed in alternate 
medium-sensitivity epochs (Table 1). Keplerian 
curves (see Figure 10) were drawn through the 
data of each epoch (circled) for ease of display. 35 
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East Offset (mas) 

Fig. 14. Map of timc-avcragcd continuum emis- 
sion obtained from a uniformly weighted average 
of images from fifteen epochs. The rms noise level 
is 0.05 mjy. Contours arc in 2.5(t steps, starting 
at 2.5(7. An average restoring beam is shown in 
the lower right. The origin of coordinates is at the 
position of the maser reference feature. 
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Fig. 15. — Light curve of fliix density for epochs 

with peak continuum emission > Scr (see Table 6). 
Strictly speaking, this is not a total flux density, 
but rather an angle integrated flux density, since 
we measure it with an interferometer that has a 
limited range of baseline lengths. 
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Fig. 16. — Position-velocity traces for red-shifted 
emission that demonstrate the magnitude of the 
postion offset affecting some epochs and the effec- 
tiveness of correction (§3.2). Both panels show the 
weighted average of red-shifted masers in medium- 
sensitivity epochs (blue) and the trace for one 
high-sensitivity observation, BM112H (red), (left) 
- No position shift to the BM112H data is ap- 
plied. The masers appear to be offset in position 
from the weighted average, (right) - Good agree- 
ment after a single position shift (applied to the 
BM112H data) that minimizes the offset in a least 
squares sense. 
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